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Transport magnetic currents driven by moving kink crystal in chiral helimagnets 
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We show that the bulk transport magnetic current is generated by the moving magnetic kink 
crystal (chiral soliton lattice) formed in the chiral helimagnet under the static magnetic field applied 
perpendicular to the helical axis. The current is caused by the non-equilibrium transport momentum 
with the kink mass being determined by the spin fluctuations around the kink crystal state. An 
emergence of the transport magnetic currents is then a consequence of the dynamical off-diagonal 
long range order along the helical axis. We derive an explicit formula for the inertial mass of the 
kink crystal and the current in the weak field limit. 

PACS numbers: Valid PACS appear here 
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How to create, transport, and manipulate spin cur- 
rents is a central problem in the multidisciplinary field 
of spintronics. 1 The key theoretical concepts there in- 
clude the current-driven spin-transfer torque^ and re- 
sultant force acting on a domain wall (DW) 3 in metal- 
lic ferromagnetic/nonmagnetic multilayers, the dissipa- 
tionless spin currents in paramagnetic spin-orbit coupled 
systems^ and magnon transport in textured magnetic 
structures. 5 A fundamental query behind the issue is how 
to describe transport magnetic currents^ Conventionally, 
the charge current is defined by the product of the carrier 
density and the drift velocity related via the continuity 
equation. In the case of spin current, the deviation of the 
spin projection from its equilibrium value plays a role of 
a charge. Then, an emergence of the transport magnetic 
currents may be expected in non-equilibrium state as a 
manifestation of the dynamical off-diagonal long range 
order (ODLRO). 7 Historically, Doring first pointed out 
that the longituidal component of the slanted magnetic 
moment inside the Bloch DW emerges as a consequence 
of translational motion of the DW4 An additional mag- 
netic energy associated with the resultant demagnetiza- 
tion field is interpreted as the kinetic energy of the wall. 

Recent progress of material synthesis sheds new light 
on this problem. In a series of magnets belonging to 
chiral space group without any rotoinversion symme- 
try elements, the crystallographic chirality gives rise to 
the asymmetric Dzyaloshinskii interaction that stabi- 
lizes either left-handed or right-handed chiral magnetic 
structures. 9 In these chiral helimagnets, magnetic field 
applied perpendicular to the helical axis stabilizes a pe- 
riodic array of DWs with definite spin chirality forming 
kink crystal or chiral soliton latticed In this paper, we 
demonstrate that the magnetic transport analogous to 
Doring effect^ occurs in the moving kink crystal of chi- 
ral helimagnets and serves an example of the dynamical 
ODLRO in non-equilibrium state. An essential point is 
that the kink crystal state has a degeneracy originat- 
ing from the translational symmetry. Consequently, the 
transport momentum has a form MX, where M and X 
represent the kink mass and the collective coordinate of 
the kink in the laboratory frame. The kink crystal be- 



haves as a heavy object with the inertial mass M. 

We start with a spin Hamiltonian describing the chiral 
helimagnet, 

H = -J Si " S i + D " J2 Sj x Sj — H • ^2 Si, 
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where the first term represents the ferromagnetic cou- 
pling with the strength J > between the nearest neigh- 
bor Si — 5* (cos Oi, sin 9% cos ipi, sin 9i sin ipi) and Sj, where 
9i and ipi denote the local polar coordinates. The sec- 
ond term represents the parity-violating Dzyaloshinskii 
interaction restricted to the nearest neighbor pairs of the 
adjacent ferromagnetic planes, characterized by the the 
mono-axial vector D = De x along a certain crystallo- 
graphic chiral axis (taken as the rc-axis). The third term 
represents the Zeeman coupling with the magnetic field 
H = 2fiBHe y applied perpendicular to the chiral axis. 
When H = 0, the long-period incommensurate helimag- 
netic structure is stabilized with the definite chirality 
(left-handed or right-handed) fixed by the direction of 
the mono- axial D- vector. 

In the continuum limit, the Hamiltonian density cor- 
responding to the lattice Hamiltonian ([I]) is written as 



H 



X - sin 2 9 (d x pf 



qo sin 2 9 (d x (p) — (3 sin 9 cos ip, 



(2) 



where the energy is measured by JS 2 , and j3 = H/JS. 
The semi-classical spin variable is represented as S = 
S (cos 9, sin cos (p, sin 9 sin ip) by using the slowly vary- 
ing polar angles 9{x) and <p(x)[see Fig. QJa)]. The 
helical pitch for the zero field (j3 = 0) is given by 
go = tan _1 (D/J) w D/J. Under the transverse field, 
a regular array of the magnetic kink is formed^ The 
kink corresponds to the phase winding in the left-handed 
(Aip = +27r) or right-handed(A<^ = — 2ir) manners. 
Since we assume the uniform mono-axial Dzyaloshinskii 
vector D, the kink with only positive (left-handed) or 
negative (right-handed) charge are energetically favored. 
The kink with the same charge repel each other, just 



like in the case of the Coulomb repulsion. So, the mag- 
netic kink crystal (soliton lattice) is formed, as shown in 
Figs. Gib) and (c). 




FIG. 1: (a) Polar coordinates in the laboratory frame, (b) 
Formation of the magnetic kink crystal in the chiral helimag- 
nets under the transverse magnetic field, and (c) concomitant 
phase modulation. In (b), we depict a linear array of the spins 
along one chiral axis that is ferromagnetically coupled to the 
neighboring arrays. 

The magnetic kink crystal phase is described by the 
stationary soliton solution, 9 = it/2 and cos (<po(x)/2) = 
sh(itix/k), where m = \/j3 corresponds to the first 
breather mass and "sn" denotes a Jacobian elliptic 
function^ The period of the kink crystal is given by 
Iq = 2kK(k)/ y/]3. The elliptic modulus k (0 < n < 
1) is determined by the energy minimization condition 
E(k)/k = irq /Am. Here, K(k) and E(k) denote the el- 
liptic integrals of the first and second kind, respectively. 

Now, we consider the fluctuations around the classi- 
cal solution and write 9(x) = it/2 + u(x), and ip(x) — 
<Po{x) + v(x). When we consider only the tangential ip- 
mode, our problem is reduced to the one first investigated 
by Sutherland?^ The ip-mode is fully studied in the con- 
text of the chiral helimagneti 13 ' 14 In the present work, 
however, it is essential to take into account not only the 
(/3-mode but the #-mode to argue the longitudinal mag- 
netic current. Expanding ([2]) up to u 2 and v 2 , we have 
H = J dx(Ho + ri u + H v + Hint)) where Ho gives the 
classical solution and H u = uC u u, H v — vC v v, where 
the differential operators are defined by 

£« = - ^(d x ip Q ) 2 + q (d x f ) + ^/3cos<p Q , (3) 



C v = --dl + i/3cos</? - 



(4) 



The lowest-order coupling between the u and v modes 
comes from Hint = —u 2 (d x v) 2 /2, that is neglected here. 
In the case of zero-field, (3 = 0, we have C u = —d 2 /2 + 
go/2, and C v — —d 2 /2. Therefore we see that u-mode 
acquires the mass go (scaled by JS 2 ), while the v-mode 
becomes massless. This situation naturally arises, be- 
cause the i;-mode is a Goldstone mode, but the u-mode 
is not. Even after switching the perpendicular field, the 
it-mode (f-niode) remains to be massive (massless). 

From now on, we argue that the massive ^-fluctuations 
carry the magnetic current. First, we perform the 
mode expansions, v{x,t) = J2 n r l n (^) Vn ( x ) an( ^ u ( x it) — 
J2 n £n(t)u n (x) , and seek the energy dispersions for 
the normal vibrational modes, satisfying C u u n (x) — 
A n u n (x), and C v v n (x) — p n v n (x), respectively. Intro- 




-2K/n 2K/n 



FIG. 2: The energy dispersions of the eigen modes for (a) the 
tangential (^-fluctuation (ujq = ^fp) and (b) the longitudinal 
^-fluctuations (ujq — VX). 



ducing x — mx/n, we have the Schrodinger-type equa- 
tions, 

d 2 u n (x)/dx 2 = [2« 2 sn 2 x 

-k 2 (1 + 2\ n /l3) - 4 + AKT]u n {x), (5) 
d 2 v n (x)/dx 2 = [2K 2 sn 2 i - k 2 (1 + 2p n //3)]v n (x), (6) 

with t = qo/m. In Eq. ([5]) we consider the case of 
weak field corresponding to small k, leading to dna: w 
1. Now, the equations ([5]) and ((6]) reduce to the 
Jacobi form of the Lame equation^ d 2 A a (x)/dx 2 = 
{£(£+ l)K 2 sn 2 x + A} A a (x), with I = 1. It is known 
that the solution is parameterized by a single continuous 
complex parameter a aS) 12 i 13 



V4 (-KX/2K) 



(7) 



where TV is a normalizing factor and di (i — 1,2,3,4) 
denote the Theta functions^ with Q being the Floquet 
indexi 12 ! 13 ' 14 The energy dispersion is obtained by deter- 
mining A = — k 2 (1+A) as a function of the Floquet index 
Q which marks eigenstates instead of n. It is known 12 
that the dispersion consists of two (generally £+1) bands 
specified by the acoustic branch A\ = k' 2 / k 2 sn 2 (a, n'), 
Qi = Tra/2KK' + Z(a,n'), and the optical branch 
A 2 = l/[K 2 sn 2 {a,K% Q 2 = TTU/2KK' + Z (a, «') + 
dn (a, k') cn (a, k') /sn (a, k'), where a e [—K',K']. 
Here, K' denotes the elliptic integral of the first kind 
with the complementary modulus k' = yl — k 2 and Z 
denotes the Zeta-function^ The complex parameter £o 
in Eq.([7|) are given by ia + K and ia for the acoustic 
and optical branches, respectively. We have the acoustic 
branch, < A 1 < k' 2 /k 2 , for < \Qi\ < tt/2K, and 
the optical branch, 1/k 2 < A 2 < oo, for tt/2K < \Q 2 \- 
The energy gap A = 1, opens at \Q\ = %/2K. We 
present the dispersions ujq in Fig. [21 where the gapless 

acoustic V psHSAy 2 and optical ^/JIbTTSA^ 2 bands 
of ^-excitations are depicted together with the gap- 
full acoustic ^/JlbHS[Ai + At/k — 4/k 2 ] 1 / 2 and optical 
\fp,BHS[A 2 + At/k - A/k 2 ] 1 / 2 bands of 6>-excitations. 

Next we consider the collective dynamics of the kink 
crystal. For this purpose, we carry out the canonical 
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formulation by using the collective coordinate method^ 
We start with the corresponding Lagrangian, 



L = J dx [5(1 



cos9) ( p t -V(9,cp)}, 



(8) 



where the Berry phase term is taken into account. The 
mode expansion leads to the vibrational term, V [9, <p\ = 
E n ^nin + J2 n PnVn- Elevating the position of the kink 
center X to a dynamical variable, we write the solution 
in the form 



<p = ifa [x - X(t)\ + Yln=l Vn(t)v n (x - X(t)) , 
9 = n/2 + Zn=iUt)u n (x-X(t)). 



(9) 



Plugging these expressions into the Lagrangian ([5]), we 
have L = SJ2 n f, n (t)K ln + SXJ2 n K 2n £ n (t) ~ 

£En,m^3nm£n(*)??m(i) ~ En A «£n ~ UnPnVl, 

with the coefficients K\ n = J dxv n (x) , K2 n 
j dx (dip /dx) u n (x), and K 3nm = J dxu m (x) v n (x) . 
This Lagrangian is singular because the determinant 
of the matrix of second derivatives of the Lagrangian 
with respect to velocities (Hessian) turns out to be 
zero. Therefore we need to construct the Hamil- 
tonian by using the Dirac's prescription for the 
constrained Hamiltonian systems. The canonical mo- 
menta conjugated to the coordinates X, £„, and r) n , i.e., 
Pi = dL/dX = SJ2 n K 2n tn, Pin = dL/d£ n = 0, p 3n = 
dL/dr] n = -SKi n -SJ2 m ^3mn£m> lead to the extended 
Hamiltonian, H* = piX + J2 n P^4n + J2 n P3nVn - L 
with a set of primary constraints, 



^ =Pi-S^2 n K 2n U = 0, 

4n = P2n=0, 
*3n = P3n + SK ln + 5 £ m #3mn£m = 0. 



(10) 



Because of a lack of primary expressible velocities the 
Hamiltonian with the imposed constraints 

2T« = ^X + Y: *£& + E ^ + H ph , 



coincides with H* , where H ph = E„ + E„/ « 7 7n- 
It governs the equations of motion of the constrained 
system, i.e. the constraints are hold at all times. This 
leads to a set of dynamical equations, 



En K 2n£,n = 0, 

-2A„£„ + XSK 2n ~ S J2 m Kznmhm = 0, 
-2p n r) n + 5 J2 m K 3mn £,m = 0, 



(11) 



to give £„ 



(2) 

constraints $„ 



and rj n = 0. Imposing the secondary 
Tj n = to be constant in time, we 
obtain f) n = 0. Together with the second constraint in 
(fTTjl this yields — (SK2 n /2X n )X , and we reach the 
final form of the physical Hamiltonian, H p h = p\/2M, 
where p\ = MX involves the soliton mass 



M 



* 2 E 



K 2 
2A n 



(12) 



Now we are ready to define the longitudinal spin cur- 
rent. We start with the linear momentum carried by the 
kink crystal, 



P = S (1- cos 9)ip x dx, 



(13) 



where Lq is the system size. By using 9 = ir/2 + u and 
tp = (fo, for a steady current, we obtain 

P « S[cp (L )-<p (0)] + S [ L °u(x)^dx 



dx 



2TrQS + Sj2tnK 2n . 



(14) 



We here introduced the topological charge, Q — [tpo(L) — 
(y9 (0)]/27r. Using the result = (SK2 n /2X n )X and 
Eq. (|12|) , we obtain an important formula, 



P = 2ttSQ + MX, 



(15) 



that plays an essential role in this paper. The first term 
is associated with the equilibrium background momen- 
tum and the second one corresponds to the transport 
current carried by the ^-fluctuations. Apparently, the 
transverse magnetic field increases a period of the kink 
crystal lattice and diminishes the topological charge Q 
and therefore it affects only the background linear mo- 
mentum. The physical momentum related with a mass 
transport due to the excitations around the kink crystal 
state is generated by the steady movement. 

The "superfluid mass current" is accompanied by the 
"superfluid magnetic current" transfered by the 9- 
fluctuations. It is determined through the definition of 
the magnetic density^ N — 5(1 — cos6>). By using 9 = 
tt/2 + u(x,t), we have dN/dt = Ssin9d9/dt w Sdu/dt, 
with u(x,t) = J^n £n(t)u n [x — X(t)\ . Therefore, for a 
steady current, we obtain the continuity equation 



-x- 



,d_ 

dx 



E 



S 2 K 2 
2X n 



dj x 
dx ' 



(16) 



where we introduced the magnon time-even current car- 
ried by the ^-fluctuations 



(17) 



Here we used the fact £ n = because of the constraint. 
The time evenness is manifested by appearance of not X 
but X 2 . The important point to note is that the only 
massive 9-mode can carry the longitudinal magnon cur- 
rent as a manifestation of ordering in non- equilibrium 
state, i.e., dynamical off-diagonal long range order. 

The final stage is to explicitly compute (jTTj) . We can 
exactly prove that the optical branch does not contribute 
to the magnetic current because of the orthogonality 
(the proof is detailed in a later paper). After a lengthy 
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but straightforward manipulation, the contribution of the 
acoustic branch is obtained as a function of x — mx/n, 



fix) 



8E 2 (k)X< 



J^ql (4E(k)/k - 1) 



dn (i) 



2X 2 



dn (i) , 
(18) 



for the weak field case corresponding to small k leading 
to E(k) ss 7r/2. On the other hand, the background spin 
currentii is shown to become jbg{x) ~ dipo(x) / dx — qo cx 
dn (i) — 2E(k)/it. We stress that the physical meaning 
of jbg is completely different from the current described 
by Eq. (|18| in non- equilibrium state.— We present a 
schematic view of an instant distribution of spins in the 
current-carrying state in Fig. [3ja). In Fig. [3jb), we 
present a snapshot of the position dependence of the cur- 
rent density j x (x) = iac(^) m the weak field limit, given 
by Eq. fiB. 



(a) 
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(b) 
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and the effect is not hindered by dephasing. Finally, we 
propose possible experimental methods to trigger off the 
spin current considered here. Spin torque mechanism: 
the spin-polarized electric current can exert torque to 
ferromagnetic moments through direct transfer of spin 
angular momentum, This effect, related with Aharonov- 
Stern effect^ for a classical motion of magnetic moment in 
an inhomogeneous magnetic field, is eligible to excite the 
sliding motion of the kink crystal by injecting the spin- 
polarized current (polarized electron beam) in the direc- 
tion cither perpendicular or oblique to the chiral axis. 
The spin current transported by the soliton lattice may 
amplify the spin current of the injected carriers. XMCD: 
to detect the longitudinal magnetic currents accompa- 
nied by the dynamical ODLRO, x-ray magnetic circular 
dichroism (XMCD) may be used. Photon angular mo- 
mentum may be aligned either parallel or anti-parallel 
to the direction of the longitudinal net magnetization. 
Ultrasound: further possibility to control and detect the 
spin current is using a coupling between spins and chiral 
torsion. 19 Ultrasound with the wavelength being adjusted 
to the period of the kink crystal may excite the periodic 
chiral torsion and resonantly supply the kinetic energy to 
the kink crystal. Consequently, the ultrasound attenua- 
tion may occur.— TOP technique: the most direct way of 
detecting the traveling magnon density may be winding 
a sample by a pick-up coil and performing the time-of- 
flight (TOF) experiment. Then, the coil should detect a 
periodic signal induced by the magnetic current. 



FIG. 3: (a) A schematic view of an instant distribution of 
spins in the current-carrying state. This picture corresponds 
to the case of intermediate field strength, (b) A snapshot of 
the position dependence of the current density j x (x) in the 
weak field limit, exactly treated in this paper. j x (x) is scaled 
by its maximum j„ lax = j x (0). 

In realizing the bulk magnetic current proposed here, 
a single crystal of chiral magnets serves as spintronics 
device. The mechanism involves no spin-orbit coupling 
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